Few-Body Systems 0, 1-6 (2008) 



Few- 
Body 
vstems 




© by Springer- Vcrlag 2008 
Printed in Austria 



Explicit Representations for the T-Matrix on 
Unphysical Energy Sheets and Resonances 
in Two- and Three-Body Systems* 



Alexander K. Motovilov 

Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 
Joliot-Curie 6, 141980 Dubna, Moscow Region, Russia 



Abstract. We describe basic structure of the two- and three-body T-matrices, 
scattering matrices, and resolvents continued to the unphysical energy sheets. 
The description is based on the explicit representations that have been found 
for analytically continued kernels of the T-operators. 



1 Introduction 

Resonance is one of the most interesting and intriguing phenomena in quantum 
scattering. With a resonance one usually associates an unstable state that only 
exists during a certain time. The original idea of interpreting resonances in quan- 
tum mechanics as complex poles of the scattering amplitude (and hence, as those 
of the scattering matrix) goes back to G. Gamov (1928). For radially symmetric 
potentials, the interpretation of two-body resonances as poles of the analytic 
continuation of the scattering matrix has been entirely elaborated in terms of 
the Jost functions. Beginning with E. C. Titchmarsh (1946) it was also realized 
that the s-matrix resonances may show up as poles of the analytic continuation 
of the Green functions. 

Another, somewhat distinct approach to resonances is known as the complex 
scaling (or complex rotation) method. The complex scaling makes it possible 
to rotate the continuous spectrum of the A^-body Hamiltonian in such a way 
that resonances in certain sectors of the complex energy plane turn into usual 
eigenvalues of the scaled Hamiltonian. In physics literature the origins of such 
an approach are traced back at least to C.Lovelace (1964). A rigorous approval 
of the complex scaling method has been done by E. Balslev and J. M. Combes 
(1971). A link between the s-matrix interpretation of resonances and its complex 
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rotation counterpart was established by G.A. Hagedorn (1979) who has proven 
that for a wide class of potentials the scaling resonances are also the scattering 
matrix resonances. 

If support of the interaction is compact, the resonances of a two-body system 
can be treated within the approach created by P. Lax and R. Phillips (1967). An 
advantage of the Lax-Phillips approach is in the opportunity of giving an elegant 
operator interpretation of resonances. The resonances show up as the discrete 
spectrum of a dissipative operator which is the generator of the compressed 
evolution semigroup. 

For more details on the history of the subject and other approaches to reso- 
nances, as well as for a bibliography we refer to paper Here we only notice 
that, in contrast to the "normal" bound and scattering states, the resonant ones 
are still a mysterious object and many questions related to resonances remain 
unanswered. In particular, it is still unknown how to describe (in an indisputable 
way) the scattering of a particle on a resonant state of two other particles. To say 
the truth, there is a problem even with definition of resonance since resonances 
are not a unitary invariant of a self-adjoint (Hermitian) operator. Moreover, 
following to B.Simon [2], one may conclude that no satisfactory definition can 
rely on a single Hamiltonian and always an extra structure is necessary. Say, 
an unperturbed dynamics (in quantum scattering theory) or geometric setup (in 
acoustical or optical problems). Resonances are always relative as the scattering 
matrix itself. 

In our approach we follow the typical setup where the resonances arising 
due to an interaction V are considered relative to the unperturbed dynamics 
described by the kinetic energy operator Hq. The resolvent G{z) = {H — z)~^ 
of the total Hamiltonian H = Hq + F is an analytic operator-valued function of 
z € C \ cr{H). The spectrum cr(H) of -ff is a natural boundary for holomorphy 
domain of G{z) considered as an operator-valued function. However the kernel 
G{-, •, z) may admit analytic continuation through the continuous spectrum of H. 
Or the form {G{z)(p, ip) may do this for any ip,ip oi a, dense subset of the Hilbert 
space Sj. Or the "compressed" resolvent PG{z)P admits such a continuation for 
P the orthogonal projection onto a subspace of ^. In any of these cases one deals 
with the Riemann surface of an analytical function. 

In the simplest example with H = Hq = —A, the two-body kinetic energy 
operator in coordinate representation, we have 

G[x,x ,z) = - j — , 
An \x — x'\ 

where x,x' are three-dimensional vectors. Clearly, G{x, function of 

the energy z has a two-sheeted Riemann surface which simply coincides with 
that of the function z^^'^. 

In this way one also arrives with the concept of the unphysical energy sheet (s). 
The copy of the complex energy plane where the resolvent G{z) is considered 
initially as an operator- valued function is called the physical sheet. The remain- 
der of the Riemann surface is assumed to consist of the unphysical sheets (in 
general, an unphysical sheet may only be a small part of the complex plane). 
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Analytic function is uniquely defined by its values on an infinite set in C 
having limiting point(s). Thus, if one knows the resolvent (T-matrix, S'-matrix) 
on the physical sheet then one may, in principle, to express it on unphysical sheets 
through its values in the physical sheet. In this case all the study of resonances 
would reduce to a work completely on the physical sheet. 

In we have found just such expressions. More precisely, we have derived 
explicit representations for the two- and three-body G{z), T{z), and S{z) on 
unphysical energy sheets in terms of these quantities themselves only taken on 
the physical sheet. In particular, the representations obtained show which blocks 
of the scattering matrix are "responsible" for resonances on a certain unphysical 
sheet. 



2 Two-body problem 

In general, we assume that the interaction potential v falls off in coordinate 
space not slower than exponentially. When studying resonances of a two-body 
system with such an interaction one can employ equally well both coordinate and 
momentum representations. However in the three-body case it is much easier for 
us to work in the momentum space (for an explanation see p. 149). This is 
one of the reasons why we proceed in the same way in the two-body case. Thus, 
for the two-body kinetic energy operator /iq we set {hof){k) = k'^f[k) where 
G M'^ stands for the reduced relative momentum. In case of a local potential 
we have v{k,k') = v{k — k') and v{k) = v{—k). For simplicity we assume that 
the function v{k) is holomorphic in k on the whole three-dimensional complex 
space C^. 

The transition operator (t-matrix) reads 

t{z) = V - vg{z)v, (1) 

where g{z) = (h — z)~^ denotes the resolvent of the perturbed Hamiltonian 
h = ho + v. The operator t is the solution of the Lippmann-Schwinger equation 

t{z) = v-vgo{z)t{z), (2) 

that is, 

t{k, k', z) = v{k, k') - [ d^ -(k,q)tig,k',z) 

taking into account that gQ{k, k' , z) = 5{k — k')/{k'^ — z). 

Clearly, all dependence of t on z is determined by the integral term on the 
r.h.s. part of © that looks like a particular case of the Cauchy type integral 

^{z)= [ dq ^ (4) 
' ' J^N ^ X + q^ - z 

for = 3. Cauchy integrals of the same form but for both = 3 and A^ = 6 we 
will also have below in three-body equations of Sec. 3. 
Let 9^A) A G C, be the Riemann surface of the function 



C(^) 



{z-Xy/^, N odd, 
log(2: — A), A^ even. 



4 Explicit Representations for the T-Matrix on Unphysical Sheets 

If N is odd, is formed of two sheets of the complex plane. One of them, where 
{z — A)^/^ coincides with the arithmetic square root \/ z — X, we denote by Uq. 
The other one, where {z — A)^/^ = —\/z — A, is denoted by iJi. 

If is even, the number of sheets of 9\x is infinite. In this case as the index 
£ of a sheet IIi we take the branch number of the function log(z — A) picked up 
from the representation log(z — A) = log l-z — A| + i 2Tri + i4> with (p E [0, 27r). 

The following statement can be easily proven by applying the residue theorem 
(if necessary, see '3^ for a proof). 

Lemma 1. For a holomorphic f{q), q € C^, the function ^{z) given by Q) is 
holomorphic on C\ [A, +oo) and admits the analytic continuation onto ^\ as 
follows 

^{z)\n, = <P{z) -£7ri{^/^^ [ dqf{VI^q), (5) 

where S^~^ denotes the unit sphere in centered at the origin. (Position of 
the argument z in the sheet Uq on the r.h.s. part of © is the same as that of 
the argument of'P\ni on the sheet Ui.) 

Now set {go{z)fi,f2) = /jg3 dq ^^i'^}'^'' where /i and /2 are holomorphic. 
Then by Lemma 1 

(5o(^;)/i,/2)|77i = (5'o(^)/i,/2)|i7o --n-iy/z / dq fi{^/zq)f2{y^q). 

which means that the continuation of the free Green function go{z) onto the 
unphysical sheet 77i can be written in short form as 

9o{z)\ni = go{z) + ao{z)j\z)j{z), (6) 

where ao(-z) = —ni^/z and j{z) is the operator forcing a (holomorphic) function 
/ to set onto the energy shell, i.e. [j{z)f){k) = f{^/zk). 

Taking into account ©, on the unphysical sheet iJi the Lippmann-Schwinger 
equation (j^J turns into 

t' = v - v{go + aoj^j)t', t' = t\ni . 

Hence {I+vg^jt' = v—a^j'^j t' . Invert I+vgo by using the fact that t{z) = v—vg^t 
and, hence, (/ + vgQ)~^v = t: 

t' = t- aotj^jt'. (7) 

Apply j{z) to both parts of Q and obtain jt' = jt — aojtj^jt', which means 

{I + ao jtj^) jt' =jt. (8) 

Observe that I + aojtj'^ is nothing but the scattering matrix s{z) since the kernel 
of s{z) reads 

s(fc, k\ z) = (5(fc, k') — t{^k., yfzk' , z). 

Hence jt' = [s{z)]^^it. Now go back to Q and get t' = t — aQ tj"''[s(z)]~^jt, that 
is, 

t{z)\n, = t{z) - ao(z) t{z)]\z)[s{z)\-^ ]{z)t{z). (9) 
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All entries on the r.h.s. part of © are on the physical sheet. This is just the 
representation for the t-matrix on the unphysical sheet we looked for. 

From one immediately derives representations for the continued resolvent, 

g{z)\n, = 5 + ao (I - gv)]^ [s{z)]~^j{I - vg), (10) 

and continued scattering matrix, 

s{z)\n, =£[s{z)]-'£, (11) 

where £ is the inversion, {£f){k) = f{—k). Hence, the resonances are nothing 
but zeros of the scattering matrix s{z) in the physical sheet. That is, the energy 
z on the unphysical sheet TTi is a resonance if and only if there is a non-zero 
vector A of L2{S'^) such that s{z)A = for the same z on the physical sheet. 

The function A{k) is the breakup amplitude of the resonance state. This 
means that in coordinate space the corresponding "Gamov vector" (the reso- 
nance solution to the Schrodinger equation) has the following asymptotics 

g-iv^|a;| 

V'res(a;) ~ A{-x) . 

x-^oo \X\ 

3 Three-body problem 

Let Hq be the three-body kinetic energy operator in the center-of-mass system. 
Assume for simplicity that there are no three-body forces and thus the total 
interaction reads V = vi + V2 + where Va, a = 1,2, 3, are the corresponding 
two-body potentials having just the same properties as in the previous section. 

The best way to proceed in the three-body case is to work with the Faddeev 
components [1] 

Ma(3 = 5apVa- VaG{z)vi3 («,/? = 1,2,3) 

of the T-operator T{z) = V — VG{z)V where, as usually, G{z) denotes the 
resolvent of the total Hamiltonian H = Hq + V. The components satisfy 
the Faddeev equations 

Mal3{z) = 5c,f3ta{z) - ta{z)Go{z) M^p{z) (12) 

with Go{z) = {Hq - 2:)"^ and ta{P, P' , z) = ta{ka, k'^, z - p'^)6{pa - p'a) where 
f^a,Pa denote the corresponding reduced Jacobi momenta (see [1^ for the precise 
definition we use) and P = {ka,pa) G M.^ is the total momentum. 

Assume that any of the two-body subsystems has only one bound state with 
the corresponding energy 6^ <i 0, a — 1,2,3. Assume in addition that all of 
these three binding energies are different. It is easy to see that the thresholds 
£i) ^2, ^3, and are associated with particular Cauchy type integrals in the 
integral equations ()12|). By Lemma 1 the two-body thresholds Ea appear to be 
square-root branching points while the three-body threshold is the logarithmic 
one. In order to enumerate the unphysical sheets we introduce the multi-index 
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I = (4,^1,^2,^3) with £0 = •••,-1,0, 1,... and £a = 0,1 if a = 1,2,3. Clearly, 
only encircling the two-body thresholds one arrives at seven unphysical sheets. 
The three-body threshold generates infinitely many unphysical sheets. (There 
might also be additional branching points on the unphysical sheets, in particular 
due to two-body resonances.) 

It turns out that the analytically continued equations H12|) can be explicitly (!) 
solved in terms of the matrix M = {Mq,^} itself taken only on the physical sheet, 
just like in the case of the two-body t-matrix. But, of course, now the result 
depends on the unphysical sheet 7T£ concerned. More precisely, the resulting 
representation reads as follows 

M\n,=M + QMLSj^LQM. (13) 

In the particular case we deal with, L and L are 4x4 scalar matrices of the 
form L = diag(£o,^i, 4,4) and L = diag(|4|,^i,4,4), respectively; St{z) = 
/ -|- L{S{z) — I)L is a truncation of the total scattering matrix S{z) and the 
entries Q m , Qm are explicitly written in terms of the half-on-shell kernels of M 
(see formula (7.34) of 1 ). From H13|) one also derives explicit representations for 
G{z)\ni and S{z)\ni, similar to those of (|in|) and (fTT|) . respectively. 

Thus, to find resonances on the sheet 11 i one should simply look for the zeros 
of the truncated scattering matrix Si{z), that is, for the points z in the physical 
sheet where equation Si{z)A = has a non-trivial solution A. The vector A will 
consist of breakup amplitudes of the resonance state into the various possible 
channels. Within such an approach one can also find the virtual states. 

In order to find the amplitudes involved in Si, one may employ any suitable 
method, for example the one of Refs. El El based on the Faddeev differential 
equations. In these works the approach we discuss has been successfully applied 
to several three-body systems. In particular, the mechanism of emerging the 
Efimov states in the ^He trimer has been studied [Hlini^ 
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